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We present a first and extensive study of thresliold resummation effects for supersymmetric 
(SUSY) particle production at hadron colliders, focusing on Drell-Yan like slepton-pair and slepton- 
sneutrino associated production. After confirming the known next-to-leading order (NLO) QCD 
corrections and generalizing the NLO SUSY-QCD corrections to the case of mixing squarks in the 
. . . virtual loop contributions, we employ the usual Mellin A^-space resummation formalism with the 

' minimal prescription for the inverse Mellin-transform and improve it by resumming 1/A-suppressed 

' and a class of A-independent universal contributions. Numerically, our results increase the theoret- 

, ical cross sections by 5 to 15% with respect to the NLO predictions and stabilize them by reducing 

■ the scale dependence from up to 20% at NLO to less than 10% with threshold resummation. 

c : 
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^ ■ I. INTRODUCTION 

■ 

I> ', The Minimal Supersymmetric Standard Model (MSSM) [ll, I3| is one of the most promising extensions of the 
I. Standard Model (SM) of particle physics. It postulates a symmetry between fermionic and bosonic degrees of freedom 

■ in nature, predicting thus the existence of a fermionic (bosonic) supersymmetric (SUSY) partner for each bosonic 
,. (fermionic) SM particle. Its main advantages are the stabilization of the gap between the Planck and the electroweak 

scale [1| , gauge coupling unification at high energy scales Q , and a stable lightest supersymmetric particle as a dark 
matter candidate [5|. Spin partners of the SM particles have not yet been observed, and in order to remain a viable 
solution to the hierarchy problem, SUSY must be broken at low energy via soft mass terms in the Lagrangian. As a 
consequence, the SUSY particles must be massive in comparison to their SM counterparts, and the Tevatron and the 
LHC will perform a conclusive search covering a wide range of masses up to the TeV scale. 
I Scalar leptons are among the lightest supersymmetric particles in many SUSY-breaking scenarios ^ and often decay 

Q-i into the corresponding Standard Model (SM) partner and the lightest stable SUSY particle. A possible signal for 
^ ' slepton-pair production at hadron colliders would thus consist in a highly energetic lepton pair and associated missing 
energy. An accurate calculation of the transverse-momentum spectrum [7j allows us to use the Cambridge (s)transverse 
^ ' mass to measure the slepton masses Q and s pin Q and to distinguish this signal from the SM background, which 
is mainly due to WW and ti production [13, UH- Current experimental (lower) limits on electron, muon, and tau 
^ slepton masses are 73 GeV, 94 GeV, and 81.9 GeV, respectivel y f]3. The leading-order (LO) cross section for the 
' production of non- mixing slepton-pairs has been calculated in jl3l.ll4ill5l[l^. while the mixing between the interaction 
eigenstates was included in ^17|]. The next-to-leading order (NLO) QCD corrections have been calculated in [l^, and 
the full SUSY-QCD corrections have been added in [l^. However, the presence of massive non-mixing squarks and 
gluinos in the loops makes the genuine SUSY corrections considerably smaller than the standard QCD ones. 

In this paper, we extend this last work by including mixing effects relevant for the squarks appearing in the 
loops, and we consider the threshold-enhanced contributions, due to soft-gluon emission from the initial state. These 
contributions arise when the initial partons have just enough energy to produce the slepton pair in the final state. 
In this case, the mismatch between virtual corrections and phase-space suppressed real-gluon emission leads to the 
appearance of large logarithmic terms a"[ln^"~^(l — z)l{l — z)]+ at the n*'^ order of perturbation theory, where 
z = M"^ js, M is the slepton-pair invariant mass, and s is the partonic center-of-mass energy. When s is close to 
Af^, the large logarithms have to be resummed, i.e. taken into account to all orders in as- The convolution of the 
partonic cross section with the steeply falling parton distributions enhances the threshold contributions even if the 
hadronic threshold is far from being reached, i.e. r — A-P / S ^ 1, where S is the hadronic center-of-mass energy. 
Large corrections are thus expected for the Drell-Yan production of a slepton pair with invariant mass M of a few 
100 GeV at the Tevatron and LHC. 
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All-order resummation is achieved through the exponentiation of the soft-gluon radiation, which does not take place 
in z-space directly, but in McUin A^-space, where N is the Mellin- variable conjugate to z and the threshold region z ^ 1 
corresponds to the limit N ^ oo. Thus, a final inverse Mellin-transform is needed in order to obtain a resumnied 
cross section in z-space. Threshold resummation for the Drell-Yan process was first performed in [2l| at the 
leading-logarithmic (LL) and next-to-leading-logarithmic (NLL) levels, corresponding to terms of the form a" In^" N 
and a" In^""^ N. The extension to the NNLL level (a" In "^^ N terms) has been carried out both for the Drell-Yan 
process [l^ and for Higgs-boson production [1^. Very recently, even the NNNLL contributions (a" ln^"~^ iV terms) 
became available [24, 25, 26]. 

In this paper, we will perform resummation for slepton-pair production at the NLL level, the reason being that 
away from the threshold region, the resummed calculation has to be matched to the fixed-order calculation, which 
is at present only known to NLO accuracy. Our analytical results for neutral (7, Z°) and charged {W^) current 
slepton-pair and slepton-sneutrino associated production at NLO will be presented in Sec. [TTl In Sec. IIIIl our fixed- 
order results will be used to perform the threshold-resummation to NLL accuracy, and numerical predictions will be 
made in Sec. IIVI We summarize our results in Sec. |Vl The mixing of sfermion interaction eigenstates is discussed in 
App. [Al and the NLO SUSY-QCD form factors for mixing squark loop contributions are collected in App. |b1 



II. SLEPTON-PAIR PRODUCTION AT FIXED ORDER IN PERTURBATIVE QCD 

In this section, we present the leading and next-to-leading order contributions in the strong coupling constant 
as to the mass-spectrum for slepton-pair and slepton-sneutrino associated production in hadronic collisions through 
neutral- and charged-current Drell-Yan type processes, 

ha{Pa)hb{pb) ^ ~k{pi)lf*{p2). (1) 

We define the square of the weak coupling constant = e^/sin^ 9w in terms of the electromagnetic fine structure 
constant a — / (47r) and the squared sine of the electroweak mixing angle xw = sin^ 9w- The coupling strengths of 
left- and right-handed (s)fermions to the neutral and charged electroweak currents are then given by 

{Lzff^Rzff} = 2Tj-2efXw, (2) 

{L^^j^^R^j^j^ ^ {LzffS{,Sf;,RzffS{,sf;}, (3) 

^wff = V2 cos OwVff, (4) 

^vv/Jj = ^wff S(j^S^^*, (5) 

where the weak isospin quantum numbers are — ±1/2 for left-handed and — for right-handed (s)fermions, 
their fractional electromagnetic charges are denoted by ey, and Vff are the usual CKM-matrix elements. In general 
SUSY-breaking models, the sfermion interaction eigenstates are not identical to the respective mass eigenstates, and 
mixing effects must be included in the coupling strengths through the unitary matrices S-^ diagonalizing the sfermion 
mass matrices (see App. |A|) . For purely left-handed sneutrino eigenstates i^l a diagonalizing matrix is not needed, i.e. 

and otherwise. In non-minimal flavour violating (NMFV) models with inter-generational mixing, the 

squark and slepton (sneutrino) mass matrices become in principle six- (three-) dimensional, although flavour-changing 
neutral lepton currents are experimentally strongly constrained [27[. 

Thanks to the QCD factorization theorem, the unpolarized hadronic cross section 
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a,b 



dXa dXbfa/h^{Xa,^J■F) fb/htiXb,^J■F)^ab[z,M^■,asi^ij{),—,— ] (6) 



can be written as the convolution of the relevant partonic cross section cTab with the universal distribution functions 
fa.b/hc,b partons a, 6 inside the hadrons /la.fc, which depend on the longitudinal momentum fractions of the two 
partons Xa^b and on the unphysical factorization scale /i^. The partonic scattering cross section will be expressed in 
terms of the SUSY particle masses rrij , mo, rriq., rrig and the masses of the neutral and charged electroweak gauge 
bosons mz and mw- The dependence on the strong coupling constant a^, the factorization and renormalization 
scales and /ijj, the invariant mass of the slepton pair M and the scaling variable z = M'^ / s, where s = XaX^S and 
S = {pa +pb)'^, will be explicitly shown. In QCD perturbation theory, the partonic cross section can be expanded in 
powers of as, 

z, M2; as[ii\), = ^ ^1;' z, M^; ■ (7) 
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FIG. 1: Feynman diagram for slepton-pair {V — 7, and slepton-sneutrino associated {V = W^) production at leading 
order in perturbative QCD. 



In this work we have computed the LO {n ~ 0) and NLO (n = 1) coefficients in the case of general mixing between 
the sfermion interaction eigenstates. 



A. Total Cross Sections at Leading Order 



At leading order in perturbative QCD, slepton-pair and associated slepton-sneutrino production proceed through 
an s-channel exchange of a photon, a Z"- or a W^^-boson, 



qq' ^W^^hvt, !*Di, 



(8) 



as shown in Fig. [TJ The corresponding cross sections were first calculated in [l^, [13, ITsl . [iGj , and mixing effects 
between the interaction eigenstates relevant for third-generation sleptons were added in [17[. 
For neutral currents, the first coefficient of Eq. ^ is given by 
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(9) 



The three terms in Eq. ([5]) represent the squared photon-contribution, the photon-Z° interference and the squared 
^^-contribution, respectively. The slepton-mass dependence is factorized in the velocity 



The purely left-handed charged-current cross section is easily derived from Eq. ^ by setting 



mz mw, 



ei — Rzq 
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(10) 



(11) 



which gives 



(0) , 2. 



a'oiAnSil-z) 
9M2 



Lwqq'L^fl^-^ 



i2xl,{\-xw?{l-ml,/M^f 



(12) 



After integration of the differential cross sections presented in [T^ [T3| , we find agreement with the neutral-current 
result of [Tsj and the charged-current result of [iGj in the limit of non-mixing mass-degenerate sleptons, obtained by 
setting all mixing matrices to the identity and by summing over the left- and right-handed eigenstates, as well as with 
the general charged- and neutral-current results of [l7| . 
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FIG. 2: Contributions of virtual diagrams for slepton-pair {V = 'y,Z'^) and slepton-sneutrino associated {V = W^) produc- 
tion at next-to- leading order in perturbative QCD. The first and second lines show the QCD and SUSY-QCD corrections, 
respectively. In the SUSY-QCD case, one has to sum over squark mass-eigenstates k, I — 1,2. 




FIG. 3: Contributions from real gluon emission diagrams for slepton-pair {V = 7, Z'^) and slepton-sneutrino associated {V — 
W^) production at next-to-leading order in perturbative QCD. 



B. Next-to-Leading Order SUSY-QCD Corrections with Mixing Squark Loops 

The NLO QCD and SUSY-QCD corrections to the slepton-pair production cross section have been studied for 
non-mixing sleptons in [Tsl . [l^ . At NLO in perturbative QCD, the quark-antiquark annihilation process receives 
contributions from virtual gluon exchange (see upper part of Fig. [2]) and real gluon emission (see Fig. [3]) diagrams, 
and we also have to take into account the quark-gluon initiated subprocess (see Fig. S]) . The infrared and coUinear 
singularities of the three-parton cross sections are extracted using the dipole subtraction formalism [2^ . and the 
virtual corrections have been evaluated in the MS renormalization scheme. For the SM QCD diagrams one has the 
well-known results uM, 12911 
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which expose the LO cross-sections crQ\AP) in factorized form. Cp — 4/3 and = 1/2 are the usual QCD cofour 
factors, and Pqq^qg are the AltareUi-Parisi sphtting functions [sO] 



Cf 
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■S{l-z) + 



(1 + ^) 



and 



(15) 
(16) 



We remind the reader that our normahzation corresponds to a perturbative expansion in powers of a^/Tr. 

The three lower diagrams of Fig. [5] contain SUSY-QCD corrections. Generalizing results from [l^ to the case of 
mixed squark mass eigenstates fc, Z = 1, 2 in the virtual loop diagrams, we obtain for neutral and charged currents 
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^(1 — z) and 
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(17) 



(18) 



where now only the diagonal squared photon contribution to the Born cross section factorizes. The virtual loop 
coefficients fj, f-yz^ fz and fw are given in App.|B] Note that the quark mass, which appears in the off-diagonal mass 
matrix elements of the squarks running in the loops, corresponds to a linear Yukawa coupling in the superpotential 
and can not be neglected, even if it is much smaller than the total center-of-mass energy of the colliding partons 
allowing for a massless factorization inside the outer hadrons. 

The full NLO contributions to the cross section are then given by 
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III. THRESHOLD RESUMMATION AT NEXT-TO-LEADING LOGARITHMIC ORDER 

In this section, we recall some well-known results about soft-gluon resummation for the Drell-Yan process. The 
hadronic cross section in Eq. ^ can be written in factorized form in McUin A^-space as 

a{N, M2) = f-/hAN + l,P%) h/h,[N + l,/i|^) aab[N, as^M^fil, M^ti^), (21) 

ab 
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where the A'^-moments of the various quantities are defined according to the MelUn transform 



(22) 



with y ~ T, z, and Xa,b, respectively, for F ^ a, a, and fa.b 



A. Soft-Gluon Resummation in Mellin A'^-Space 



The terms leading to finite and singular contributions in the threshold region are those proportional to 5(1 — z) and 
to ln(l — z) and the plus-distributions in Eqs. and (O. They respectively give rise to constant (A^-independent) 
and to (lnA^)/A^ and In' A^(i — 1,2) terms in Mellin space. [For a derivation of these correspondences we refer the 
reader, for instance, to App. A of [1^.] These are the contributions that need to be resummed to all orders in 
perturbative QCD. In the following, we suppress the dependence on M^//i|j and M^/mf for brevity. 

The moments of the partonic cross section can be written in resummed form as [20, l2l| 



^T'^\N,as) = a^" Cabias) exp 5(iV,a,) 



The A'^-independent terms are collected in the Cab functions 

oc 



(23) 



(24) 



(25) 



These contributions are mostly due to hard virtual corrections to the cross section, i.e. the terms proportional to 
(5(1 — 2) in Eqs. and p^ . The exponential form factor S can be written at NLL as 



S{N,as)^2 / dz— / 

Jo ./mI 



(26) 



The integrand in Eq. (|26p embodies the contributions coming from the coUinear emission of soft gluons from initial- 
state partons, i.e. the terms proportional to the plus-distributions in Eq. (|13p . It is a series expansion in the strong 
coupling constant, 



00 



(27) 



whose coefficients are perturbatively computable through a fixed-order calculation. In particular, it has been proven 
[3l| that in the MS factorization scheme the coefficients of the A-function are exactly equal to the large- coefficients 
of the diagonal splitting function 



lqq{as) = dzZ 
JQ 



Pgg{z) = -A{a,) \nN + 0{l), 



(28) 



where N = N exp[7E] and 75=0.5772... is the Euler-Mascheroni constant. Performing the integration in Eq. 
and using Eq. ((27| . we obtain the form factor up to NLL, 

SiN,as)^giiX) lnAr + 52(A). 

The functions gi and 172 resum the LL (a" In"^^ A^) and NLL (a" In" A^) contributions, respectively, and are given by 

MM 



(29) 



gi(A) = ^ [2A+(l-2A)ln(l-2A)] and 
PoA 



52 (A) 



PI 



2A + ln(l - 2 A) + 2 (1 - 2 A) 



A2 

PI 



[2A + ln(l-2A)] 



+ ^[2A + ln(l-2A)]ln-^-— ^In-^, 
Po I^R Po 



(30) 



(31) 
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where A = [Pocts \nN]/7r. The first two coefficients of the QCD /3-function are 

Po^^{nCA-2Nf) and pi ^ ^{17 Cl - 5 Ca Nf ~ 3Cf Nf), 

Nf being the number of effectively massless quark flavours and Cp = 4/3, Ca = 3 the usual QCD colour factors. 
Thus, the knowledge of the first two coefficients of the function A{as) [33. l33l|. 



Ai = Cf and A2 = ^Cp 



Ca 



67 tt' 



18 6 J 9 

together with the first coefficients of the C-functions in Eqs. ([M]) and (P5|) . 



Nf 



allows us to perform resummation up to NLL. 
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and 



(32) 

i. 

(33) 

(34) 
(35) 



B. Improvements of the Resummation Formalism 

In the limit of large TV, the cross section is clearly dominated by terms of 0(ln^ N), ^^(In A^) and 0{1). It seems thus 
reasonable to neglect terms suppressed by powers of 1/N in the resummation formalism. Actually these last terms 
are multiplied by powers of In and could as well provide a non- negligible effect in the threshold limit. In [ssl. i36] it 
has been shown that these contributions are due to coUinear parton emission and can be consistently included in the 
resummation formula, leading to a "coUinear-improved" resummation formalism. The modification simply amounts 



to the introduction of an iV-dependent term in the C'^Jw and Cqg^ coefficient of Eqs. and ([55)1 
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2 fi-i 
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(36) 
(37) 



Furthermore, the exponentiation of the contributions embodied in the C-function has been proved in [37[, leading 
to the following modification in Eq. (p3|) : 



a 



(RES) 
ab 



exp 



exp 



SiN,as 



(38) 



As the authors of Ref. [33| recognize, this exponentiation of the A^- independent terms is not comparable to the standard 
threshold resummation in terms of predictive power. While in the latter case a low-order calculation can be used to 
predict the behaviour of full towers of logarithms, in the former case it is not possible to directly get information on 
the behaviour of constant terms at, say, n loops, but a complete calculation at the n**^ perturbative order will still be 
necessary. Nonetheless, the comparison of the numerical results obtained with and without the exponentiation of the 
constant terms can at least provide an estimate of the errors due to missing higher-order corrections. 



C. Inverse Mellin- Transform and Matching Procedure 



Once resummation has been achieved in TV-space, an inverse Mellin-transform back to the physical cc-space is 
needed. The customary way to perform this inversion, avoiding the singularities of the A^-moments, is the "Minimal 
Prescription" of [33 |. 



1 fCMP+ioo /A/f2\ 

' dNi'^) a{N,M^). 



a = 



2"^* J Cm p-i 00 



S 



-N 



(39) 



The constant Cmp has to be chosen so that all the poles in the integrand are to the left of the integration contour 
in the complex A^-plane except for the Landau pole at A^ = exp[7r/(2 (3q as)], which should lie far to the right on the 
real axis. 
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Finally, a matching procedure of the NLL resummed cross section to the NLO result has to be performed in order 
to keep the full information contained in the fixed-order calculation and to avoid possible double-counting of the 
logarithmic enhanced contributions. A correct matching is achieved through 



.(F.O.; 



1 



2711 



CMP+i oo 



dN 



M2 



(40) 



where a^^^'^■^ is the fixed-order perturbative result, cr^^^^) is the resummed cross section, and cr^^-^^' is the truncation 
of the resummed cross section to the same perturbative order as a^^ '^'K In our case the expansion of the resummed 
partonic cross section up to order as reads 



^(--)(iV,M^) = aW(Af^) 1 + ^ C.(21n^^-21niVln^)+C;;i 

In Mellin-space, the fixed order NLO cross sections of Eqs. (fT3|) and (fT4|) read [sH] 
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+ 0{al) and 

TV^ + 11 iV3 22 + 14 TV + 4 
iV2 {N + lf{N + 2f 



N{N + 1){N + 2) /4 



with Si{N) = '£f=i In the large-TV limit, Si{N) ~ IniV 1/(2 ^^), and we get 
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2 In^iV- 2 In TV In 
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+ c% 



0[al) and 



TT «f 



(41) 
(42) 



(43) 



(44) 



(45) 
(46) 



Comparing Eqs. ([41]) . (|42|) . (|45)) . and (|46p . we see that the expansion of the resummed cross section at order as 
correctly reproduces the fixed order result in the large- A'^ limit, including even terms that are suppressed by 1/7V. 



IV. NUMERICAL RESULTS 



For the masses and widths of the electroweak gauge bosons, we use the current values of mz = 91.1876 GeV, 
mw = 80.403 GeV, Tz = 2.4952 GeV, and Tw = 2.141 GeV. The CKM-matrix elements are computed using 

(C12C13 S12C13 siae""* \ 

-.S12C23 - Ci2S23Sl3e"' C12C23 - 3125235136"^ S23C13 (47) 
S12S23 - Ci2C23Sl3e*'^ -C12S23 - Si2C23Sl3e*'^ C23C13 / 

with Sij — sin Oij and = cos Oij , 9ij being the usual angles relative to the mixing of two specific generations i and 
i and 5 being the CP-violating complex phase. Their average values are given by 

S12 = 0.2243 , S23 = 0.0413 , S13 = 0.0037 , and (5 = 1.05. (48) 

The squared sine of the electroweak mixing angle 

sin^6'vy = l-m^/m% (49) 
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and the electromagnetic fine structure constant 

a ^ ^PlGFrri^wSm^Ow/T: (50) 

can be calculated in the improved Born approximation using the world average value oi Gp — 1.16637 • 10^^ GeV"'^ 
for Fermi's coupling constant ^12]. 

The physical masses of the SUSY particles and the mixing angles are computed with the computer program SUS- 
PECT [3^, including a consistent calculation of the Higgs mass, with all one-loop and the dominant two- loop radiative 
corrections in the renormalization group equations that link the restricted set of SUSY-breaking parameters at the 
gauge coupling unification scale to the complete set of observable SUSY masses and mixing angles at the electroweak 
scale. We choose one minimal supergravity (mSUGRA) point, SPS la, and one gauge-mediated supersymmetry 
breaking (GMSB) point, SPS 7, as benchmarks for our numerical study SPS la is a typical mSUGRA point 

with an intermediate value of tan/3 = 10 and > 0. It has a model line attached to it, which is specified by 
Too = — j4o = 0.4 TOi/2- For TOi/2 = 250 GeV, this SUSY-breaking scenario leads to light sleptons fi, ei, f2, 62, 
Vt and i>e with masses of 136.2, 146.4, 216.3, 212.3, 196.1 and 197.1 GeV and to heavy squarks with masses around 
500-600 GeV [i^. SPS 7 is a GMSB scenario with a fi as the next-to-lightest SUSY particle (NLSP) and an effective 
SUSY-breaking scale A = 40 TeV, A^mos = 3 messenger fields of mass Mmes = 80 TeV, tan/3 = 15, and ^ > 0, which 
leads again to light sleptons with masses of 114.8, 121.1, 263.9, 262.1, 249.5 and 249.9 GeV, respectively, and even 
heavier squarks with masses around 800-900 GeV. Its model line is defined by Mmes 2 A. The slepton masses and 
mixing angles are actually quite similar for the SPS la mSUGRA and SPS 7 GMSB points, so that the corresponding 
production cross sections will not differ significantly. Slepton detection will, however, be slightly different in both 
scenarios, as the sleptons decay to a relatively massive neutralino (x^) lightest SUSY particle (LSP) at SPS la, but 
to a very light gravitino LSP at SPS 7. The lightest tau slepton thus decays into a tau lepton and missing transverse 
energy. Feasibility studies of tau-slepton identification at the LHC with the ATLAS detector [4l[ and tau tagging 
with the CMS detector have recently shown that stau masses should be observable up to the TeV range. 

Our cross sections are calculated for the Tevatron pp-collider, currently operating at ^/S ~ 1.96 TeV, as well as for 
the LHC pp-coUider, bound to operate at ^fS = 14 TeV starting in 2008. For the LO (NLO and NLL) predictions, 
we use the LO 2001 jill (NLO 2004 ^) MRST-sets of parton distribution functions (PDFs). For the NLO and 
NLL predictions, as is evaluated with the corresponding value of A^^^ = 255 MeV at two-loop accuracy. We fix the 
unphysical scales ^.p and equal to the invariant mass M of the slepton (slepton-sneutrino) pair. 

A. Invariant-Mass Distributions for Slepton Pairs 

The invariant-mass distribution M^Aa/AAI for first- (and equal-mass second-) generation sleptons at the Tevatron 
is shown in Fig. [5] (top), the one for (slightly lighter) third-generation sleptons at the LHC in Fig. [5] (bottom). In both 
cases, we have chosen the SPS 7 GMSB benchmark point. The differential cross section dcr/dM has been multiplied 
by a factor AI^ in order to remove the leading mass dependence of propagator and phase space factors. As is to 
be expected for P-wave production of scalar particles, the distributions rise above the threshold at y/S = 2to; with 
the third power of the slepton velocity /3, see Eq. (|9]), and peak at about 100 GeV above threshold (at 370 GeV for 
M^da/dM and 310 GeV for da /AM for the Tevatron; 410 GeV and 300 GeV for the LHC), before falling off steeply 
due to the s-channel propagator and the decreasing parton luminosity. As can also be seen from Eqs. ([T^ and p^ . 
the QCD corrections do not alter the P-wave velocity dependence close to threshold. At the Tevatron, the total and 
NLO SUSY-QCD predictions exceed the maximal LO cross section by 36 and 31%, respectively, whereas at the LHC, 
the maximal cross section increases by 28 and 27%. Threshold resummation effects are thus clearly more important 
at the Tevatron, where the hadronic center-of-mass energy is limited and the scaling variable t = M"^ / S is closer to 
one, and they increase with M to the right of both plots. The maximal theoretical error is estimated in Fig. [5] by 
an independent variation of the factorization and renormalization scales between M/2 and 2M. It is indicated as a 
shaded, vertically, and horizontally hatched band for the LO, NLO SUSY-QCD, and the total prediction. At LO, the 
only dependence comes from the factorization scale. It increases with the momentum-fraction x of the partons in the 
proton or anti-proton and is therefore already substantial for small M at the Tevatron, but only for larger M at the 
LHC. At NLO, this dependence is reduced due to the factorization of initial-state singularities, but a strong additional 
dependence is introduced by the renormalization scale in the coupling q;s(^_r). After resummation, this dependence 
is reduced as well, so that the total scale uncertainty at the Tevatron diminishes from 20%-35% for NLO to only 
16%-17% for the matched resummed result. The reduction is, of course, more important in the large-M region. At 
the LHC, where as is evaluated at a larger renormalization scale and is thus less sensitive to it, the corresponding 
numbers are 18%-25% and 15%-17%. 
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FIG. 5: Invariant-mass distribution M^da/dM of e^-pairs at the Tevatron (top) and fi-pairs at the LHC (bottom) for the 
benchmark point SPS 7. We show the total NLL+NLO matched and the fixed order NLO SUSY-QCD and LO QCD results, 
including the respective scale uncertainties as horizontally hatched, vertically hatched and shaded bands. 
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FIG. 6: isT-factors as defined in Eq. (|51|l for et-pair (top) and associated ez,!/* production (bottom) at the Tevatron for the 
benchmark point SPS la. We show the total NLL+NLO matched result, which is almost identical to the purely resummed 
result at NLL, as well as the fixed order NLO SUSY-QCD and QCD results. The latter practically coincides with the resummed 
result expanded up to NLO. 
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FIG. 7: A'-factors as defined in Eq. (ISlfl for fi-pair (top) and associated fii>r production (bottom) at tlie LHC for tlie benciimark 
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For the mSUGRA scenario SPS la, with nonetheless similar slepton masses and mixing angles (see above), we show 
in Figs. [5] and [7] the cross section correction factors 



where i labels the corrections induced by NLO QCD (Eqs. and additional NLO SUSY-QCD (Eqs. ^ 

and (PO)) ). resummation (Eqs. (^5]) . ([5^ . and ([57|) ). and the matched total contributions (Eq. (201)) as well as the 
fixed-order expansion (Eqs. ((4T|) and ([42])) of the resummation contribution as a function of the invariant mass M. 
As one can see immediately, the mass-dependence of these corrections for charged-current associated production of 
sleptons and sneutrinos (lower parts of Figs. [S]and Fig. [7]) does not differ substantially from the mass-dependence of 
the neutral-current production of slepton-pairs (upper parts) . 

At the Tevatron (Fig. [S]), where we are close to the threshold, resummation effects are already important at low 
M (4%) and increase to sizeable 16% at large M. The NLO QCD result is thus dominated by large logarithms 
and coincides with the expanded result at the permille level. In addition, the relative importance of the (finite) 
SUSY-QCD contributions is reduced, and the total prediction coincides with the resummed prediction, since fixed- 
order and expanded contributions cancel each other in Eq. (1401) . We have also verified that exponentiating the finite 

(iV- independent) terms collected in the coefficient function C'^J),) , as proposed in [33| (see Eq. [55)) . leads only to a 
0.6%-0.8% increase of the matched resummed result. The Tevatron being a pp-collider, the total cross section is 
dominated by gg-annihilation, and (/^-scattering contributes at most 1% at small M (or small x), where the gluon 
density is still appreciable. Integration over M leads to total cross sections for the neutral (charged) current processes 
in Fig. El of 4.12 (3.92) fb in LO, 5.3 (4.96) fb in NLO (SUSY-)QCD, and 5.55 (5.28) fb for the matched resummed 
calculation. The corresponding (global) X-factors 



a' _ / dM da' /dM 

JdMdcrLO/dM 



are then 1.29 (1.27) at fixed-order and 1.35 (1.35) with resummation. 

At the LHC (Fig. [7]), sleptons can be produced with relatively small invariant mass M compared to the total 
available center-of-mass energy a/S", so that z = rjixaXh) = M'^ / s <C 1 and the resummation of (1 — z)-logarithms is 
less important. This is particularly true for the production of the light mass-eigenstates of mixing third-generation 
sleptons, as shown in Fig. [71 In the low-M (left) parts of these plots, the total result is less than 0.5% larger than the 
NLO (SUSY-)QCD result. Only at large M the logarithms become important and lead to a 7% increase of the if-factor 
with resummation over the fixed-order result. In this region, the resummed result approaches the total prediction, 
since the NLO QCD calculation is dominated by large logarithms and approaches the expanded resummed result. 
However, we are still far from the hadronic threshold region, so that both resummed and fixed-order contributions 
and a consistent matching of the two are needed. At low M, where finite terms dominate, the resummed contribution 
is close to its fixed-order expansion and disappears with M . In the intermediate-M region, one can observe the effect 
of SUSY-QCD contributions, in particular the one coming from the qqg-vertex correction (lower left diagram in Fig. 
[5]). As M > 2mg, one crosses the threshold for squark-pair production and observes a resonance in Fig. [7] As for the 

Tevatron, exponentiating the finite (iV-independent) terms collected in the coefficient function C'^^m leads only to a 
1% increase of the matched resummed result. The LHC being a high-energy pp-collider, it has a significant gluon- 
luminosity, in particular at small M (or x), and indeed the gg-subprocess changes (lowers) the total cross section by 
7% at small M and 3% at large M. After integration over M, we obtain total cross sections of 27 (9.59) fb in LO, 
34.3 (12.3) fb in NLO SUSY-QCD, and 34.6 (12.5) fb for the resummed-improved result, corresponding to global 
ii'-factors of 1.28 for fixed-order and 1.29 for the matched resummed cross section for both processes. Resummation 
of large logarithms is thus not as important as for the Tevatron at the benchmark point SPS la. 



B. Scale Variations of the Total Cross Section 



In this section, we study the dependence of total, i.e. invariant-mass integrated, slepton-pair and slepton-sneutrino 
cross sections on three different scales: first the dependence on the effective SUSY-breaking scale A as defined in the 
GMSB model line SPS 7, and then the dependence on the renormalization and factorization scales fJ-R,F at the point 
SPS 7. We remind the reader that the model line for SPS 7 is M^cs — "2 A with Ny^cs — 3, tan/3 — 15, and /x > 
fixed. The benchmark point at A = 40 TeV will be indicated in the figures where A varies as a vertical dashed line, 
and we show in addition to the scale A the mass scale of the produced charged slepton, ranging from 80 (87.5) to 280 
(385) GeV for (n). 
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In Fig. [H total cross sections (top) and i^- factors (bottom) are shown at the Tevatron, which is expected to produce 
a total integrated luminosity of 4-8 fb~^. Only left-handed (charged) eigenstates couple to the weak (electromagnetic) 
neutral current, so that their respective cross sections are enhanced. However, even for e^-pair production the mass- 
range is limited to masses below 280 GeV, where the cross section reaches 0.1 fb and at most one event would be 
produced. NLO and resummation corrections are clearly important, as they increase the LO prediction by 18 to 28% 
(lower part of Fig. [8]). At the SPS 7 benchmark point, the corrections would thus induce a shift in the selectron mass 
as deduced from a total cross section measurement by about 8 GeV (cf. the two dashed lines in the upper part of Fig. 
|5]). By comparing the NLO and total predictions, one observes an increased importance of threshold resummation for 
heavier sleptons, as expected. 

At the LHC (Fig. [51), even very heavy sleptons can be detected (see above), but we restrict ourselves to the range 
A < 135 TeV. The NLO and resummed corrections are again large (25-30%), but the resummation corrections only 
become appreciable for large SUSY- breaking scales (or slepton masses). The largest cross section is obtained for pair 
production of the light stau mass eigenstate, even though it has a large right-handed component. Conversely, the 
heavier stau mass eigenstate has a large left-handed component, so that its cross section is less suppressed. At the 
SPS 7 benchmark point, the corrections would again induce a shift in the slepton (fi) mass as deduced from a total 
cross section measurement by about 8 GeV (cf. the two dashed lines in the upper part of Fig.[5|). 

Finally, we consider the theoretical uncertainty of invariant-mass integrated total cross sections at the Tevatron (Fig. 
[TU]) and the LHC (Fig. [TT]) as induced by variations of the factorization scale (top), renormalization scale (middle), 
or both (bottom). The /i^-dependence (middle), which is absent in LO, is first introduced in NLO, but then tamed 
by the resummation procedure. On the other hand, the logarithmic /ii?-dependence (top), already present through 
the PDFs at LO, is overcompensated (reduced) at NLO for the LHC (Tevatron) and then (further) stabilized by 
resummation. This works considerably better at the LHC, where at least one quark PDF is sea-like and the PDFs are 
evaluated at lower x, than at the Tevatron, where both PDFs can be valence-like and are evaluated at relatively large 
X. In total, the theoretical uncertainty at the Tevatron (LHC), defined by the ratio of the cross section difference at 
f^F = IJ-B. = 'nT'il'^ and jip = A*_r = 2m^ over their sum, increases from 20 (7) % in LO to 29 (17) % in NLO, but is 
then reduced again to 23 (8) % for the resummed-improved prediction. 



In summary, we have presented a first and extensive study on threshold resummation effects for SUSY-particle 
production at hadron colliders, focusing on Drell-Yan like slepton-pair and slepton-sneutrino associated production 
in mSUGRA and GMSB scenarios. After confirming the known NLO QCD corrections and generalizing the NLO 
SUSY-QCD corrections to the case of mixing squarks in the virtual loop contributions, we employed the usual Mellin 
iV-space resummation formalism with the minimal prescription for the inverse Mellin-transform, but improved it by 
resumming 1/A^-suppressed and a class of A'^-independent universal contributions. Numerically, our results increase 
the theoretical cross sections by 5 to 15% with respect to the NLO prediction and stabilize them by reducing the scale 
dependence from up to 20% at NLO to less than 10% with threshold resummation. 
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The soft SUSY-breaking terms Af oi the trilinear Higgs-sfermion-sfcrmion interaction and the off-diagonal Higgs 
mass parameter ^ in the MSSM Lagrangian induce mixings of the left- and right-handed sfermion eigenstates fL,R of 
the electroweak interaction into mass eigenstates /i,2- The sfermion mass matrix is given by Q 
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FIG. 9: Total cross sections (top) and _K'-factors as defined in Eq. (|52p (bottom) for tliird-generation slepton-pair and slepton- 
sneutrino associated production at the LHC along the model line attached to the SPS 7 benchmark point (vertical dashed line). 
We show the total NLL+NLO matched and the fixed order NLO (SUSY-)QCD and LO QCD results. 
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FIG. 10: Dependence of the total cross section for first- (and second-) generation slepton pairs at the Tevatron on the factor- 
ization scale (top), renormalization scale (middle), and both scales (bottom) for the SPS 7 benchmark point. 
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m 



LL = + {Tj - e/ sin^ Ow) m'z cos 2/3, (A2) 



= mp,+ef sin Ow'Tn^ cos 2/3, (A3) 

^ _ * f cot /3 for up — type sfermions. (A4) 

i^LB. f IJ- Wa,n/3 for down — type sfermions. ^ ' 

It is diagonalized by a unitary matrix , Ai^ S^"^ = diag (m^, 7122), and has the squared mass eigenvalues 

™?,2 = ™/ + ^ + T - to|j^)2 + 4to2 l^i^pj . (A5) 

For real values of rriLi^, the sfermion mixing angle ^j, < 0j < 7r/2, in 
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FIG. 11: Dependence of the total cross section for third-generation slepton pairs at the LHC on the factorization scale (top), 
renormalization scale (middle), and both scales (bottom) for the SPS 7 benchmark point. 



can be obtained from 

2 rrif niLR 



tan 20.^ -■■-/-^^ _ 



If m^i?, is complex, one may first choose a suitable phase rotation /jj = e^'^fa to make the mass matrix real and then 
diagonalize it for and tan/? = Vu/vd is the (real) ratio of the vacuum expectation values of the two Higgs 
fields, which couple to the up-type and down-type (s)fermions. The soft SUSY-breaking mass terms for left- and 
right-handed sfermions are nip and mp, respectively. 



APPENDIX B: SUSY-QCD FORM FACTORS FOR MIXED SQUARK MASS EIGENSTATES 



At next-to-leading order in perturbative QCD, SUSY-QCD form factors are induced by one-loop diagrams involving 
gluinos and (generally mixed) squark mass eigenstates. Those appearing in the neutral-current cross section of Eq. 
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and the one appearing in the charged-current cross section of Eq. (jlSp by 
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The functions Bof{p'^, m\^rn^, Bq^{p'^ ^m\, m\) and CQf{p\, {pi +^2)^,^21 "^ii "^2' "^3) finite parts of the scalar 

two- and three-point functions 
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Our results agree with those of Ref. |45i] in the case of mass-degenerate non- mixing squarks. 
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